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Abstract 

We present a scheme of analytical calculations determining the criti- 
cal temperature and the number of condensed atoms of ideal gas Bose- 
Einstein condensation in external potentials with ID, 2D or 3D periodicity. 
In particular we show that the width of the lowest energy band appears as 
the main parameter determining the critical temperature of condensation. 
Is obtained a very simple, proportional to 1/3 degree, regularity for this 
dependence. The fundamental role of tunneling in physics of condensate 
establishment is underscored. 

1 Introduction 

Degenerate Bose gases provide an excellent ground for the theoretical study of 
quantum fluids since their diluteness makes possible first-principles approaches [1] 
Thanks to today's atomic physics powerful experimental techniques their prop- 
erties can be studied quantitatively through a wide range of temperature and 
densities. In the last few years Bose-Einstein condensates (BEC) and Fermi 
gases in optical lattices have been an extremely active area of research[2]. Pe- 
riodic potentials have been used to examine the transport of Bose-condensed 
samples [3], to investigate effects correlated with the physics of strongly corre- 
lated many-body systems [4]. These systems can be possibly implemented in 
quantum information processing [5] and there are proposals on how to build 
quantum gates [6], [7] and qubit buses [8] for information exchange. 

The presence of trapping potentials has a big impact on characteristics of 
the condensate, and in particular on the value of temperature at which the gas 
passes from normal into BEC state and a macroscopic order coherence starts to 
form. The predicted increase of critical temperature in a trap played a big role 



in obtaining the BEC state with the help of laser cooling and evaporative cooling 
techniques. However presence of the external potential doesn't always bring to 
the increase of critical temperature of condensation (7b). In [9] was shown that 
if the external potential is periodic, Tc decreases and for asymptotically deep 
potentials tends to 0. The decrease of critical temperature for relatively deep 
potentials, created as off-resonant standing waves, was observed experimentally 
in [10]. We also gave a descriptive physical interpretation for the dependence 
of Tc on external trapping potential parameters [11], basing on behavior of the 
distance between the low-laying energy levels of translational atom. For strictly 
periodic potentials this role is naturally played by widths of the low energy 
bands, distances between them and most of all the width of the lowest energy 
band. 

In this paper, using some nondrastic mathematical approximations, we de- 
rive an elementary analytic expression for critical temperature Tc as a function 
of optical lattice parameters. The potential periodicity is assumed in one, two 
or three independent directions. The motion for the left directions in first two 
cases is assumed to be free. The analytical result justifies, in particular, our 
physical reasonings brought in [11] about the Tc behavior relative to lattice 
depth adiabatic changes, and concretizes the form of that behavior. The ob- 
tained results, we hope, will help in experimentally more accurately determining 
the range of parameters, where the ultracold Bose gases behave as ideal. 



2 The statistical problem of ideal Bose-gas in 
periodic potentials 

The impact of laser radiation fields on translational motion of atoms is presented 
as a momentum exchange in photon absorption and emission processes. This 
means that, in general, the atom (molecule) translational state evolution is 
connected with the evolution of internal states, leading to notion of potential 
for each energy level. In the case of large resonance detunings, however, it 
becomes possible to introduce the idea of potential for the center of mass, as 
for unstructured particles. Just such a situation will be assumed later, taking 
the laser field in form of standing waves creating the periodic potential. The 
coordinate space in our calculations is ordinary, three dimensional, while the 
periodic potential is present on one, two, or in all three directions. 

The principal grand canonical relation, relating the chemical potential [i 
with the number of atoms N in a system of volume V[12], in case of interest 
takes the form 

dP x dP Y \ dP z , p - , p - , p . = N. 

(1) 

Here Pj (j = X, Y, Z) is the atomic momentum for free motion directions and is 
the quasimomentum for the periodic potential directions. Momentum (quasi- 
momentum) is scaled by the recoil momentum 2hk, where k is the wavevector of 
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the countcrpropagating waves forming the standing wave. The energies Sj(Pj), 
as well as the chemical potential \i and the thermal energy kbT, are scaled in 
recoil energy units Er — (2fifc) 2 /2M, corresponding to momentum value 2hk. 
Formulae (1) assumes an extended form for the dispersion relation between 
energy and quasimomentum. 



2.1 ID periodic potential 



After taking the elementary integrals over two free directions in (1), we get 
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(2) 



where the integration over the total range of quasimomentum is replaced by two 
half-range integrals - from zero to infinity. To proceed, it is convenient to divide 
the whole range of integration into a sequence of Brillouin zones and expand 
the integrand logarithm into convergent Taylor series. Then the above relation 
takes the following form: 
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Now let's temporarily concentrate on the integral term and take a dispersion 
relation, for example, for a biparabolic form of the periodic potential, introduced 
in [11]: 

cos(2^P) = 1 + 2Gn(£)G 22 (e), (4) 

(the explicit expressions of Gu(e) and G 22(e) are not necessary for our later 
presentation) . 

As is seen from Fig.l, the dependence of energy e on quasimomentum P, 
calculated by relation (4), is almost linear in frame of the first energy band 
and this linearity isn't rapidly lost for higher energy bands. Such a behavior 
prompts us to introduce a new, main in context of this paper, approximation 
taking the e(P) -dependence inside each energy band as linear. Here we come 
from the fact, that in assumed thermal equilibrium state of the Bose gas the 
population and respective contribution of upper energy bands into the left-side 
value of (3) decreases quite rapidly. 

After denoting boundary energy values of the m-th zone by and 
and performing the mentioned linearization, the dispersion relation (4) will have 
the form 
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(5) 
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and the integral in (3) is expressed by zero order modified Bessel function Iq(x) 
with an exponential factor. As a result we obtain the following, more simple 
form for the main statistical relation: 
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(6) 
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where — (e^x + £ < ^il)/2 is the mean energy in the m-th band and the 
parameter = £max — in the argument of modified Bessel function is 
the width of that band. The contribution of each energy band in this relation 
is now determined by the expression in curve brackets and presents a single- 
variable convergent series, a very convenient form for numerical calculations. 
To be convinced in convergent nature of the mentioned series, one would ad- 
dress to asymptotic formulae of the modified Bessel function for great values of 
argument: 

The asymptotic behavior of the mentioned sum is then determined by expression 

1 
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(8) 



which evidently converges for any [i < ej™^ as usual for bosonic systems. 

Thus, from mathematical viewpoint relation (6) is reasonably defined and 
is convenient for calculating the functional dependence of [i on T with any 
in advance prescribed accuracy. In particular the critical temperature will be 
decided by substituting fi = e^- n : 




m—1 \n— 1 



n B T c J \ 2n B T c J 



From general principles of the Bose -Einstein condensation theory directly 
follows that for T < Tc in (6) we should interpret TV as the number of noncon- 
densed atoms(A r „ c ) and not as the total number, simultaneously taking [i = e^ ln 
which corresponds to the critical temperature. The number of condensed atoms 
is determined by the complementary relation N c — N — N nc . 

Now let's go back to formulae (7) and use it not for showing the convergence 
of the series, but for getting a new, much more simple approximate form of 
problem solution. Really, the formulae (7) can be used for small values of n 
too, if first of all n B T < 5^/2, that is if atom's thermal energy is appreciably 
smaller than the half- width of the first, most narrow energy band (for critical 
and lower temperatures this condition is well satisfied in today's experiments). 
Nevertheless, in general, the fulfillment of this condition is not mandatory and 
the use of (7) for the first addends too would be allowed approximation when 
the impact of these first members is small. Without going into details we will 
just mention that for the above mentioned replacement is enough to satisfy the 
condition 

ir 2 N 2k b T 
k 3 Vn B T > ( ' 

After substituting (7) into (6) the sum over n can be calculated and expressed 
by confluent hypergeometric function $(a, 3/2; 1). Then the main statistical 
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relation takes the form 



l,3jr 00 1 / ( m ) \ 

J-^ T ?' 2 E i exp $ (e-(^-^,3/2;l) = AT, 

m=l \ / 

(11) 

from which the critical temperature is decided as before, substituting /1 = £^ n : 



(12) 

First term in curve brackets stands for the number of atoms in the first energy 
band. The number of atoms in any, higher laying energy band is determined 
by the respective addend of the series: m = 2 for the second (first excited) one, 
and so on. 

Often in experiments with optical lattices we can limit ourselves with the 
first energy band. For this cases we get elementary relation and for making 
its appearance more contensivc, it is appropriate to introduce a notion of the 
first energy band for free (when the energy gaps tend to zero) Bose gas. If we 
denote this width by (^(in normalized units it is equal to 1/4), and the critical 
temperature of a free ideal gas [12] by Tco , than we arrive at the following final 
form: 



Ml) 



2/3 



x(l) 



T c o- (12) 



This formulae is one of the main results of the present paper. The width 
5^ is the only parameter depending on external potential. The deepening of 
the periodic potential, as is well known, compresses this width and, cosequently, 
decreases Tc, the critical temperature of condensation (such a result, based on 
numerical calculations and physical reasonings, we have presented in [11]). The 
formulae shows that this decrease happens by a very simple law: proportional 
to 1/3 degree of the first energy band width. 

The approximate character of formulae (13) is seen in fact, that in the limit 
of free gas (6^ — > 6^) the critical temperature Tc tends not to Tc but 
to (7r/2) 2 / 3 w 1.35 times higher value. For not very deep potentials (in E r 
units smaller than one) the main approximation done here is the substitution of 
dispersion curves with straight lines. The deepening of the potential (decreasing 
<5^ m ') straightens the dispersion curves, starting from the first energy band, and 
consequently softens the role of this approximation. But now grows the role 
of using the asymptotic formulae (7), as for the first members of series the 
condition /2kbT > 1 fails. The single band approximation for ultracold 

gases has a minor impact whatever the case. 



5 



2.2 2D and 3D periodic potentials 

Suppose that 2D (3D) periodic potential is a sum of two (three) periodic po- 
tentials, each one of which is periodic only in one direction. The calculation 
scheme doesn't undergo any qualitative changes and we will write the final re- 
sults straight away. Instead of (6), we come to relation 
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for 3D periodicity, where the new notations are the full analogies of the ID case. 
After applying (7) these relations are simplified and take the forms 
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correspondingly. For critical temperatures we get just as simple relations as was 
(13): 



for potentials with 2D and 
with 3D periodicity. 
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3 Conclusions 



Formulas (12), (13), (14) and (15) represent the main results of this paper. 
They show that with good approximation the width of the first energy band 
is the only parameter determining the ideal gas Bose-Einstein condensation 
critical temperature in a field of periodic potentials. From the Bloch theory 
of periodic potentials is well known that the deepening of the potential rapidly 
narrows energy bands and especially the first one. Therefore even the 1/3 
degree proportionality in ID case will be enough for decreasing the critical 
temperature twice for 2Er — 3Er potential depths. The comparison of the 
above mentioned formulas also shows that in a single-band approximation each 
direction of periodicity acts as though independent and by that speeds up the 
critical temperature dependence on the potential depth when the periodicity 
passes from ID to 2D and 3D correspondingly. 

The band structure of energy spectrum assums the existence of Bloch-type 
wavefunctions, the modulus of which is periodically spread over the whole po- 
tential. For the states which are of interest, that is with energies smaller than 
the potential height, the only mechanism of the Bloch state realization is the 
quantum tunneling through the potential barriers. By this we conclude that the 
phenomena of tunneling has an exceptional role for the obtained BEC critical 
temperature behavior regularity. If we exclude the possibility of tunneling for 
quantum particles (in our case of atom), then each low-energy atom would've 
stay trapped and localized only in one well-type region of the potential, as for 
example in harmonic potential case. The atom translational motion energy spec- 
trum would be discrete and the energy levels would move away as the potential 
depth increases. Then it would bring to an increase of the critical temperature, 
a result well known for isolated wells and just opposite to the above obtained 
ones. 
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Figure 1: The dependence of the right hand side(RHS) of (4) from normalized 
energy e. The optical lattice depth is taken IEr. The bold intervals of the 
curve correspond to allowed energy bands. 
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